Invisibility in non-Hermitian tight-binding lattices 
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Reflectionless defects in Hermitian tight-binding lattices, synthesized by the intertwining operator 
technique of supersymmetric quantum mechanics, are generally not invisible and time-of-flight mea- 
surements could reveal the existence of the defects. Here it is shown that, in a certain class of 
non-Hermitian tight-binding lattices with complex hopping amplitudes, defects in the lattice can 
appear fully invisible to an outside observer. The synthesized non-Hermitian lattices with invisible 
defects possess a real-valued energy spectrum, however they lack of parity-time (VT~) symmetry, 
which does not play any role in the present work. 

PACS numbers: 03.65.-w, 71.10.Fd, 42.82.Et, 72.20.Ee 



I. INTRODUCTION 



In recent years, the subject of invisibility physics has 
attracted a great and renewed interest, mainly triggered 
by the publication of a few seminal papers by Pendry and 
Leonhardt on transformation optics and electromagnetic 
cloaking 0, Q , which has lead to the first experimental 
observation of invisibility at microwave frequencies 
Since then, a large body of works inspired by the 
concepts of transformation optics has been published, 
and applications to matter wave cloaking have been 
suggested as well An invisible object or scatter 

is, by definition, an object which does not scatter any 
wave incident upon it; that is, a wave which shines 
on the object is not reflected or absorbed, but instead 
it is transmitted in such a way that it appears to the 
outside observer as if there were no object present. The 
concepts and methods of invisibility based on the idea of 
transformation optics apply to two- or three-dimensional 
objects. In one-dimensional systems, the possibility 
to achieve an invisible scatter is closely related to the 
realization of reflectionless potentials. For continuous 
media, this problem was investigated in a pioneering 
work by Kay and Moses in 1956 [5], and then studied 
in great detail in the context of the inverse scattering 
theory @, 0] and supersymmetric quantum mechanics 
for Hermitian systems [8|. The potentials obtained by 
such techniques, though being transparent, are generally 
not invisible. This is due to the dependence of the phase 
of the transmitted wave on energy, which is generally 
responsible for some delay and/or for the distortion of a 
wave packet transmitted across the potential [9] . 
The possibility of synthesizing reflectionless potentials 
has been also investigated for wave scattering on a 
lattice, in which wave transport occurs due to hopping 
among adjacent sites of the lattice. In the mathematical 
literature, this problem is solved by the inverse spectral 
theory of Jacobi operators, i.e., second order symmetric 
difference operators in this context, Darboux 



transformations and the intertwining operator technique 
of supersymmetric quantum mechanics have been suc- 
cessfully extended to the discrete Schrodinger equation, 
with applications to the synthesis of transparent (i.e. 
reflectionless) defects in Hermitian tight-binding lattices 
[Til Il2j . An optical realization of a special class of these 
reflectionless potentials on a lattice has been recently 
proposed for waveguide arrays and coupled-resonator 
structures with modulated coupling rates [13| . suggest- 
ing new possibilities for pulse and beam shaping. For 
Hermitian lattices, such reflectionless potentials are 
nevertheless not invisible because the bound states of 
the lattice modify the time-of-flight of a wave packet and 
generally also distorts its shape: the existence of defects 
in the lattice, though being transparent, could be then 
inferred form simple time-of-flight measurements. 
It is the aim of this work to show that fully invisibility 
of localized defects can be realized in non-Hermitian 
tight-binding lattices, which are synthesized by iterated 
application of the intertwining operator technique 
(Darboux transformation) to a defect-free tight-binding 
Hermitian lattice. The study of non-Hermitian tight- 
binding lattices has received in recent years a great 
attention (see, e.g., (T3 - [l9l | and references therein); such 
previous studies have been mainly focused to lattices 
possessing parity-time (VT) symmetry and were framed 
in the context of non-Hermitian quantum mechanics 
[17H20l |. however the possibility to realize invisibility in 
a non-Hermitian lattice was not investigated in such 
previous works [2lj |. It should be noted that the class of 
non-Hermtian lattices synthesized in the present work by 
application of the Darboux transformation and showing 
the property of invisibility are not P7~-symmetric. Nev- 
ertheless, their energy spectrum is real-valued because 
they are isospectral to an Hermitian lattice. Therefore, 
VT symmetry does not play any role in the realization 
of invisible defects discussed in this work. 
The paper is organized as follows. In Sees. II and III, 
the intertwining operator technique and its application 
to the synthesis of tight-binding lattices with reflec- 
tionless defects are briefly reviewed. The scattering 
and invisibility properties of the synthesized lattices 
are discussed in Section IV; in particular, it is shown 
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that, as for any Hermitian lattice invisibility can never 
be achieved and time-of-flight measurements can be 
used to reveal the existence of defects in the lattice, in 
non-Hermitian lattices with certain complex hopping 
rates invisibility can occur. The main conclusions are 
outlined in Sec.V, whereas some mathematical details 
and a possible realization of non-Hermitian lattice 
models based on light propagation in optical waveguide 
arrays are presented in three Appendixes. 



not belong to the point spectrum nor to the continu- 
ous spectrum of Hi. More precisely, <pn is given by an 
arbitrary superposition of two linearly-independent solu- 
tions to Eq.(2), which behave asymptotically as <pn ~ 

exp(±win) at n — > ±00 for /ii > 2k, or as (f>n ~ 
(— 1)™ exp(±win) at n — > ±00 for /zi < —2k, where 
uji > is the root of the equation 2KCOsh(wi) = 
It can be then shown by direct calculations that the fol- 
lowing factorization for Hi holds 



Hi = Qifti + Ml 



II. THE INTERTWINING OPERATOR 
TECHNIQUE FOR SPECTRAL ENGINEERING 
OF TIGHT-BINDING LATTICES 

The synthesis of reflectionless and invisible defects in 
a tight-binding lattice discussed in the next sections is 
based on the discrete analogs of the intertwining opera- 
tor technique of supersymmetric quantum mechanics Q . 
Extensions of the intertwining operator technique to the 
discrete Schrodinger equation, together with the related 
issue of inverse scattering for Jacobi operators, have been 
discussed mainly in the mathematical literature (see, for 
instance, |10l - [l2l ]). however they are not so common in 
the physical contexts. In this section we thus provide a 
brief review of the intertwining operator technique and 
its application to the problem of spectral engineering of 
tight-binding lattices. 

Let us consider a one-dimensional tight-binding lattice 
described by the Hamiltonian 



where 
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]K n (\n-l)(n\ + \n){n-l\) + J2V n \n)(n\ (1) 
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(7) 

(8) 
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Let us then introduce the new Hamiltonian H2 obtained 
from Hi by interchanging the operators TZ\ and Q\, i.e. 
let us set 



where \n) is a Wannier state localized at site n of the 
lattice, K n is the hopping rate between sites \n — 1) and 
|n), and V n is the energy of Wannier state |n). Note that 
H turns out to be Hermitian provided that the hopping 
amplitudes n n and site energies V n are real- valued param- 
eters. Let us indicate by Hi the tight-binding Hamilto- 
nian defined by Eq.(l) with hopping amplitudes and site 

energies given by «4 an d Vn X \ respectively, and let us 
assume that dn* — ► K > and Vn — > as n — > ±00, i.e. 
that the lattice is asymptotically homogeneous and free of 
defects. Let — a c U a p be the spectrum of Hi, which 
comprises the continuous spectrum a c (the tight-binding 
band —k<E<k) and the point spectrum <j p . Our goal 
is to synthesize a new tight-binding lattice Hamiltonian 
Ti.2 of the form of Eq.(l), whose spectrum is the same 
as that of Hi , except for the addition of a new real- valued 
energy level [i\ in the point spectrum, with \\x\ \ > 2k. To 

this aim, let us indicate by l^ 1 -*} = J2n \ n ) a solution 
to the second-order difference equation 

+ + WW = Ml^ (2) 



with the asymptotic behavior 



00 for n — > ±00. 



Note that such a solution does exist because \i\ does 



h 2 = n x Q x + 



(10) 



H2 will be referred to as the partner Hamiltonian of Hi. 
Using Eqs.(4-9), from E.(10) it can be readily shown that 
H2 describes the Hamiltonian of a tight-binding lattice 
[i.e., it is of the form (1)] with hopping amplitudes and 
site energies {Kn\V^} given by 



,( 2 ) = 



„(!) 



y(2) = y(l) + 
r n. r n. 



b {1) 



,(1). 



(11) 

(12) 



Note that, owing to the asymptotic behavior of «4 , Vn^ 



and (filP at n — > ±00, one has Kn' — > K and Vn' J —> for 
n — > ±00, i.e. the partner lattice described by the Hamil- 
tonian H2 is still a homogeneous lattice without defects 
at n — > ±00. An interesting property of the Hamiltonian 
H2 is that its spectrum is given by — en 1 ) U{/Ui}, 
i.e. it is the same as that of Hi plus the additional energy 
level /ii in the point spectrum. In fact, let us indicate 
by \^e) = 2n i } n(E)\ n ) a proper (or improper) eigen- 
function of Hi with energy E. Note that, if E belongs 



r(2) 
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to the point spectrum of Hi, \ip n (E)\ — > as n — > ±00, 
whereas if E belongs to the continuous spectrum of "H^ 
\ilj n (E) \ remains bounded as n — >• ±00. Since /ii does not 
belong to the point spectrum of Hi, one has E H\. Us- 
ing the factorization (3) for Hi, the eigenvalue equation 
Hi\i/je) = E\\Pe) reads explicitly 

QiftilVe) = {E - »i)\ip E ) (13) 

from which it follows that IZi^e) 7^ since E 7^ fii. 
Applying the operator IZi to both sides of Eq.(13), one 
obtains 

KiQi$b) = (E - ih)$ b }, (14) 

i.e. H2\4>e) — E\^e), where we have set \^e) = T^-i^e) 
or, explicitly [see Eq.(5)] 

ME) = r^ME) + fllHn-iiE). (15) 

Therefore, \iPe) is an eigenfunction of H2 correspond- 
ing to the energy E. Also, from Eqs.(6), (7), (15) and 
from the assumed asymptotic behavior of Kn and Vn^ 
as n — > ±00, it follows that \i))e) is a proper (improper) 
eigenfunction of H2 in the same way as V'-e) is a proper 
(improper) eigenfunction of Hi. In a similar way, one 
can show that any eigenvalue E of H2, belonging to its 
continuous or to its point spectrum, is also an eigenvalue 
of Hi provided that E ^ fn. Therefore the continuous 
and point spectra of Hi and H2 do coincide, apart from 
the energy level E — fii which needs a separate analysis. 
For E = /xi, the eigenvalue equation = Ml IV*) can 

be satisfied by taking Qi\ip) — 0, which reads explicitly 

3£tyn + ?£ l ty»+i=0. (16) 

Using the expressions of q n ^ and q n ^ given by Eqs.(6-9), 
the difference equation (16) for ip n can be solved in a 
closed form, yielding 



s/^mWi 

In view of the asymptotic behaviors of cf) n and i^ti as 
ri . — > ±00 and assuming that (j>n does not vanish for any 
integer n, it turns out that ip n is bounded and ip n — > 
as ?i — > ±00, i.e. E = jii belongs to the point spectrum 
of H2 and its eigenfunction is given by Eq.(17). 
It should be noted that the synthesis of the partner 
Hamiltonian H2, with spectrum er 2 = o"i U is not 

unique because of some freedom left in the choice of </>n^ 
satisfying Eq.(2) once /ii has been fixed: different choices 

of 4>n^ lead in fact to different lattice realizations of H2, 

(2) 

i.e. different values of hopping amplitudes K n and site 
energies Vn ■ 

The factorization method can be iterated to synthesize 
new Hamiltonians H3, Hi, H5, ... whose energy spectra 



differ from that of Hi owing to the addition of the dis- 
crete energy levels {/Ji,/i 2 }, {^1,^2,^3}, {^1,^2,^3,^4}, 

with |/i fc | >2n(k = 1,2,3,4,...). 
An interesting property, that is proven in the Appendix 
A, is the following one. Let us assume = for 

the lattice Hamiltonian Hi. Then a partner Hamilto- 
nian H2N+1, obtained from Hi by adding 2A new energy 
levels fii, fi2, fJ.3, — -,V2N with (12 = -fii, Hi = -ft,--, 
fj,2N = ~ P>2N— l) can be synthesized in such a way that 
y(2N+i) _ q This means that the partner lattice de- 
scribed by H2N+1 and supporting 2 A bound states differs 
from the original one, defined by Hi, because of different 
hopping rates K n between adjacent sites, but not for the 
site energies V n . 

As a final note, it should be mentioned that the technique 
of intertwining operators so far described could generate 
non-Hermitian lattice Hamiltonians with complex- valued 
hopping rates K n or site energies V n , even though the ini- 
tial Hamiltonian Hi is Hermitian. However, in spite of 
non-Hermiticity, the energy spectrum of such synthesized 
Hamiltonians remains by construction real- valued. This 
situation is especially interesting for the synthesis of in- 
visible defects in the lattice, as discussed in Sec. IV. 

III. TIGHT BINDING LATTICES WITH 
REFLECTIONLESS DEFECTS 

The intertwining operator technique presented in the 
previous section can be applied to the synthesis of lat- 
tices with reflectionless defects. Previous works have so 
far limited to consider Hermitian lattices (see, for in- 
stance, [ill, [HI); conversely, here we do not necessarily 
require that the partner Hamiltonians H2, H3, Hi, 
obtained by the iterated application of intertwining op- 
erator method, be self-adjoint. Notably, it will be shown 
in the next section that a truly invisibility of the de- 
fects requires the synthesis of non-Hermitian lattices. In 
this section, we first discuss the scattering properties of 
partner lattice Hamiltonians obtained by the intertwin- 
ing operator technique, and then apply the results to the 
synthesis of reflectionless defects in the lattices. 

A. Scattering properties of partner lattice 
Hamiltonians 

Let Hi and H2 be the Hamiltonians of the two part- 
ner tight-binding lattices defined by Eqs.(3) and (10). By 
construction, the two Hamiltonians have the same energy 
spectrum, except for an additional energy level fii for H2- 
The two lattices are homogeneous (i.e., free of defects) at 
n — > ±00; therefore, asymptotically they admit of plane- 
wave solutions of the form ~ exp(±«qn), where q is the 
wave number that varies in the interval < q < n. Such 
plane waves belong to the common continuous spectrum 
of the Hamiltonians, with energy E(q) = 2ncos(q). The 
reflection (ri(q), ^(q)) and transmission (ti(q), £2(9)) 
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coefficients of the two lattices are defined by the asymp- 
totic behavior of scattered waves at n — » ±00 from a 
forward-incident plane wave ~ exp(— iqn) according to 
the relations [23 



, 1 ) 1 exp(— iqn) + r 1 (q) exp(iqn) n — > —00 
ti(q) exp(-iqn) n — > 00 



for Hi, and 



ij)^J exp(-iqn) + r 2 {q) exp(ign) n 
(2(3) exp(— iqn) n 



00 



(18) 



(19) 



for "H2- Let us indicate by wi the real- valued and positive 
solution to the equation 



2k cosh(wi) 



(20) 



and let Si — (ii/\[ii\ (i.e. Si = 1 for /xj > 0, ft = —1 
for /ii < 0). It can be then proven that the following 
relations between transmission and reflection coefficients 
of the two partner Hamiltonians hold 

, / x , , s exp(-o;i/2) - ft exp(o;i/2 + iq) 

h(q) - ii(?) exp(wi/2) _ 6i exp( _ Wl/2 + iq) ( 21 ) 

exp(o;i/2) - ft exp(-u^/2 - ig) 

r 2 [q) =n{q) — j— — . (22) 

exp(u;i/2) — ft exp(— + iq) 

The proof of Eqs.(21) and (22) is given in the Appendix 
B. Here we just noticed that = |r^ 2 ^(g)| and 

I^H?) — I* (<?)l>i-e. the transmittance and reflectance 
coefficients of the two partner lattices are the same. It 
should be noted that, as |i (1) (<?)| 2 + \r m (q)\ 2 = 1 for 
the Hermitian Hi lattice, it follows that \t^(q)\ 2 + 



..(2) 



(q)\ 



1 either, even if the partner Hamiltonian 



H2 is non-Hermitian. This result is a non-trivial one 
because it is known that unitarity of the scattering ma- 
trix in a generic non-Hermitian Hamiltonian is usually 
broken, and the reflection and transmission coefficients 
can be unbounded (see, for instance, [23| and references 
therein). 

By simple iteration, Eqs.(21) and (22) can be readily 
extended to the case of the partner Hamiltonian H n ob- 
tained from Hi by adding the energy levels \ii, [X2, 
/1 at. The reflection (r]y(q)) and transmission (tjv(g)) co- 
efficients of the lattice described by Hn are given by 



, / \ , , N -rr exp(-cj fe /2) - 5k exp(cj fe /2 + iq) 
tN(q) = ti(q) M — . 

±± exp(w fc /2) - ft exp(-w fe /2 + iq) 

V exp(w jt /2) - S k cxp(-w fc /2 - iq) 



r N (q) = n(q) J] 



k=l 



exp(wfc/2) - 5k exp(-w/c/2 + iq) 



(23) 
(24) 



where Uk is the positive root of the equation 
2/icosh(w fc ) = \fi k \ and ft. = fi k /\Hk\ (k = 1, 2, 3, N). 



B. Lattice with reflectionless defects 

RefLectionless lattices containing localized defects are 
readily synthesized by assuming for Hi the Hamiltonian 



of a homogeneous and defect-free lattice (n 



(i) 



V, 



(i) _ 



0), for which ri(q) = and ti(q) = 1. In fact, 
from Eq.(24) it follows that the reflection coefficient 
rjv(g) of any partner Hamiltonian Hn vanishes, and the 
incident wave is fully transmitted through the lattice. 
Depending on the choice of the sequences 4>n \ \ <f>n , 
the resulting partner Hamiltonian may be or may 
not be Hermitian. 



Hermitian Lattices 



Examples of reflectionless and Hermitian lattices 



(a) 



(b) 
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wave number q (71 units) 
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wave number q (jt units) 



FIG. 1: (color online) (a) Behavior of the hopping rates K n 
for a Hermitian lattice as predicted by Eq.(30) for parame- 
ter values TV = 3, wi = 0.6 and a = 0. (b) Behavior of the 
hopping rates K n for a non-Hermitian lattice as predicted by 
Eq.(35) for parameter values N = 3, wi = 0.01 and a — 0.5. 
In the figure, the dots refer to Re(K n ), whereas the squares 
to Im(Kn). In (c) and (d) the behaviors of the phase of the 
transmission coefficient t(q) of the two lattices are also de- 
picted. 

obtained by the application of the intertwining operator 
technique or by other techniques have been previously 
presented in [lll - Fl3j j . The simplest case corresponds 
to the addition of a single energy level [ii outside the 
tight-binding band — k < E < n. Assuming for instance 
Hi > k, Eq.(2) can be satisfied with the choice 



4>n = cosh[wi(n — a)] 



(25) 



which ensures the Hermiticity of the partner Hamiltonian 
Hi- In Eq.(25), ui\ = acosh(/ii/2ft) and a is an arbitrary 
real parameter. The hopping amplitudes and site ener- 
gies of the partner lattice read explicitly [see Eqs.(ll) 
and (12)] 



,(2) = 



^cosh[^i(n -a-2)] cosh[^I(n - a)} 

cosh[wi(n — a — 1] 
cosh[wi(n — a + 1)] cosh[uji(n — a)] 



sh[wi(n — a)] 



sh[wi(n — a — 1)] 



(27) 
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Such a lattice, in spite of the presence of defects, is re- 
flectionless and supports one bound state, given by [see 
Eq.(17)] 



1pn = 



1 



-\/cosh[cji(n — a)] cosh[wi(n — a — 1)] 



(28) 



Another example, which was recently proposed in 
Ref . [l3j . is provided by the partner lattice H3 obtained 
from the defect-free lattice 'Hi by ad ding the couple of 
energy levels [1% > k and /12 = — £*i [24]. In this case, 
assuming again for <^>„ the expression given by Eq.(25), 
according to the analysis of Sec. II and Appendix A the 
hopping amplitudes of the Hermitian lattice H3 read ex- 
plicitly [see Eq.(A6)] 



/ cosh[wi(n — a)] cosh[<x>i(n — a — 3)] 
y cosh[wi(n — a — 1)] cosh[o;i(ri — a — 2)] 



(29) 



(3) 

whereas Vn — for the site energies. The lattice H3 is, 
by construction, rcflectionless and supports two bound 
states. With the procedure outlined in the previous sec- 
tion, Hermitian lattices supporting an arbitrarily large 
number of bound states can be constructed in this way. 
A simple and noteworthy case, which generalizes the pre- 
vious example, is provided by the lattice Hamiltonian 
H2N+1 obtained from the defect-free lattice Hi by adding 
the 2N energy levels \i\ = 2/-ccosh(u; 1 ), ^2 = —fit, ^3 = 
2k cosh(2wi), /14 = — /U3, /X2JV-1 = 2k cosh(7Vwi), 
M2JV = —^2N-l- In this case, with the choice (25) for 
<f>n , one can show that the hopping rates of the lattice 
H2N+1 take the simple form 



(2JV+1) 



cosh[o;i(n — a)] cosh[cji(n — a — 2N — 1)] 



y cosh[o;i(n — a — N)] cosh[cJi(n — a — N — 1)] 

(30) 

which generalizes Eq.(29). An an example, Fig. 1(a) 
shows the behavior of the hopping rates k n , as predicted 
by Eq.(30), for the case N — 3 and for uii = 0.6, 
a = 0. As shown in the next section, even though being 
reflectionless, such Hermitian lattices are not invisible 
owing to the energy-dependence introduced by the 
bound states in the phase of the transmission coefficient. 



N on- Hermitian Lattices 



A different choice of the sequences 
can be used to synthesized reflectionless non-Hermitian 
lattices. The simplest case corresponds, as in the pre- 
vious Hermitian case, to the addition of a single energy 
level (i± outside the tight-binding band — n < E < k. 
Let us assume, for the sake of definiteness, /ii > n and 
let us make the choice [which replaces Eq.(25)] 



= sinh[wi(n - a)], 



(31) 



where ui\ = acosh(/zi/2«;) and a is an arbitrary real (but 
non-integer) parameter. The expressions of hopping am- 



plitudes and site energies of the partner lattice Hamilto- 
nian %2 are then given by 



.(2) _ 



1/(2) = 



/sinh[wi(n — a — 2)] sinh[aji(n — a)] 



y sinh 2 [o;i(n — a — 1)] 

sinh[aji(n — a + 1)] sinh[wi(n 
sinh[wi(n — a)] 



■a)] 



ih[wi(n — a — 1)] 



(32) 
(33) 



which replace Eqs.(26) and (27), respectively. Note that, 

(2) 

as the site energies Vn are always real-valued, the hop- 

(2) . (2) 

ping amplitudes k„ are not. Specifically, n n becomes 

purely imaginary at the two lattice sites n satisfying the 
condition a < n < 2 + a. Therefore, the partner Hamil- 
tonian %2 is not Hermitian, in spite its spectrum is real- 
valued by construction. 

As a second example, let us synthesize the partner lattice 
H3 obtained from the defect-free lattice Hi by adding the 
couple of energy levels \i\ > n and ^2 = — fJ.%, assuming 
again for (f>n the expression given by Eq.(31). Accord- 
ing to the analysis of Sec. II and Appendix A, the hopping 
amplitudes of the lattice H3 now read explicitly [compare 
with Eq.(29)] 



/ sinh[aji(n — a)] sinh[wi(n — a — 3)] 
y sinh[wi(n — a — 1)] sinh[cJi(n — a — 2)] 



(3) 

whereas Vn — for the site energies. By construction, 
the lattice Hamiltonian H3 is reflectionless, has a real- 
valued energy spectrum and supports two bound states, 
corresponding to the energies E = ±2k cosh(a;i). How- 
ever, an inspection of Eq.(34) indicates that H3 is not 

(3) 

Hermitian because the hopping amplitudes Kh, take an 
imaginary value at the two sites n satisfying the condi- 
tion a < n < 1 + a and 2 + a<n<3 + a. More gen- 
erally, with the choice (31), the non-Hermitian Hamil- 
tonian H2N+1 admitting 2N bound states with ener- 
gies /ii = 2kcos1i(wi), /i 2 = — (ii, /13 = 2KC0sh(2wi), 
M4 = -A*3, .... ) H2N-1 = 2KC0sh(7Vu;i), ^ 2N = -{12N-1 
can be synthesized, corresponding to the hopping ampli- 
tudes [compare with Eq.(30)] 



(34) 



„(2N+1) _ 



/ sinh[wi(n — a)] sinh[cji(n — a — 2N — 1)] 
Y sinh[o;i(n - a - N)} sinh[o;i(n - a - N - 1)] 

(35) 

and site energies Vn N = 0. Note that the hopping 
amplitudes are purely imaginary at lattice sites n satis- 
fying the conditions a < n < a + N and a + N + 1 < 
n < a + 2N +1. As an example, Fig. 1(b) shows the 
behavior of the real and imaginary parts of the hopping 
amplitudes K n as given by Eq.(35) for N = 3, uj = 0.01 
and a — 0.5. 

One could wonder whether non-Hermitian tight-binding 
lattices with imaginary hopping amplitudes may describe 
wave transport in some physically realizable systems. 
Coupled optical waveguide structures with gain and/or 
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70. 



0.2 0.4 0.6 0.8 1 
wave number q (71 units) 

FIG. 2: (color online) Behavior of the phase ipk{q) versus the 
wave number q (in units of 71"), defined by Eq.(37), for (a) 
8k = 1, (b) <5fe = — 1 and for increasing values of ujk- The 
curves refer to the values ujk = 0.01, 0.1, 0.5, 1, 2 and 4. The 
arrows in the figures show the direction of increasing ujk- 



loss regions have been recently proposed as experimen- 
tally accessible systems to mimic the dynamics of non- 
Hermitian lattices with complex-valued site energies (see, 
for instance, [l7Hl"9l Hlj]); however, the non-Hermitian 
lattices discussed in the previous examples require imag- 
inary values of the hopping rates at some site energies, an 
issue which was not considered in such previous works. 
In the Appendix C, it is shown that suitable longitudinal 
modulations of gain/loss and propagation constants in 
evanescently-coupled optical waveguide arrays lead to an 
effective non-Hermitian lattice with imaginary hopping 
amplitudes that realizes the models discussed in this sec- 
tion. 



IV. INVISIBILITY IN NON-HERMITIAN 
LATTICES 



For a reflectionless lattice synthesized by the intertwin- 
ing operator technique, the transmission coefficient as a 
function of the wave number q of the incident wave has 
the form t(q) — exp[iLp(q)], where according to Eq.(23) 
the phase <p(q) is given by the sum of N contributions 
associated to each of the N bound states with energies 
Mi: M2, fJ-N, i-e. 



N 



<P{0) = ^2<Pk(q) 



(36) 
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FIG. 3: (color online) Propagation of an initial Gaussian- 
shaped wave packet (snapshot of the site occupation proba- 
bilities P n (t) — \{ifj(t)\n)\ 2 ) (a) in the Hermitian lattice with 
hopping rates shown in Fig. 1(a) [parameter values are given 
in the text], and (b) in the defect-free lattice. In (c) the be- 
haviors of site occupation probabilities P n (t) at time t — 70 
in the two lattices are depicted (the solid line refers to the 
Hermitian lattice with defects, the dashed line to the defect- 
free lattice). Note the advancement experienced by the wave 
packet propagating in the lattice with defects. Such an ad- 
vancement is basically ascribable to the increase of hopping 
rates k„ in the defect region [see Fig. 1(a)]. 



where 

exp[iip k (q)] 



exp(-w fc /2) - 8k exp(w fe /2 + iq) 



exp(wfc/ 2 ) - Sk exp(-wfc/2 + iq) 
2nSk cosh(cjfc), Sk = /ife/|/ifc|, and LOk > (fc 

• • 1 



(37) 




fc=i 



FIG. 4: (color online) Behavior of the hopping rates k„ for 
the Hermitian lattice, as given by Eq.(30), for iV = 3, a = 
and for decreasing values of wi (curve 1: wi = 0.6; curve 2: 
uji — 0.3; curve 3: wi = 0.2; curve 4: wi = 0.02). 

1, 2, 3, N). The behavior of <fik{q), for increasing val- 
ues of uik and for Sk = ±1, is shown in Fig. 2. In case 
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§k = 1 [i.e. /ifc > 0, see Fig. 2(a)], one has <pk (q) — 7r + q 
for uj k ;» 1 and Pk{q) —> (mod 2ir, q ^ 0) for uj k — > + . 
Similarly, in case <5fc = — 1 [i.e. Uk < 0, see Fig. 2(b)], one 
has tfk(q) — q for <^fc 1 and <Pk(q) — > (g 7^ 7r) for 
— > + . Note that, according to Eq.(36), the behavior 
of the overall phase tp(q) is given by the superposition of 
the various terms <pk (q) and does not depend on whether 
the synthesized partner Hamiltonian Hat is Hermitian or 
non-Hermitian. 

We now ask ourselves whether the defects in the partner 
lattice, in addition of being reflectionless, are also invis- 
ible to an outside observer. This condition requires that 
the phase (p(q) of the transmission coefficient be flat, i.e. 
that (dp/dq) — almost everywhere. If this condition is 
not satisfied, the spectral components of a wave packet 
crossing the defect region of the partner lattice would 
acquire the additional phase contribution ip(q), absent 
in the defect-free lattice, which would be responsible for 
a different time-of-flight and for a different distortion of 
the wave packet as compared to the same wave packet 
propagating in the ideal defect-free lattice. Therefore, 
an outside observer could detect the existence of defects 
somewhere in the lattice by e.g. simple time-of-flight 
measurements. The advance in the time of flight expe- 
rienced by the wave packet propagating in the partner 
lattice with defects can be readily calculated by stan- 
dard methods of phase or group-delay time analysis, and 
reads 



1 



dip 
dq 



1 (dip 
2k sin(<jo) V dq 



(38) 



</(i 



where qo is the carrier wave number of the wave packet 
and v g = 2Ksin(g ) > its group velocity. In particular, 
for a partner lattice synthesized by taking uj^ 3> 1, one 
has (d(p(q)/dq) ~ N (see Fig. 2), and thus the advance- 
ment of the wave packet measured by an outside observer 
(i.e. far from the defect region) would be ~ N/v g . Hence, 
comparing the time of flight measurements in the two lat- 
tices, the observer can estimate the number N of bound 
states of the partner lattice. From the above considera- 
tions, it follows that the necessary and sufficient condi- 
tion for a reflectionless lattice to be also invisible is that 
ojk — > 0. For a Hermitian lattice, from Eqs.(26) and (27) 
[and similarly from Eqs.(29) or (30)] it follows that in 
this limit the lattice becomes defect-free, i.e. K n — » 1 
and V n —¥ regardless of the value of the parameter a. 
This means that, for the Hermitian lattices synthesized 
in Sec. III. B, the invisibility condition is the absence of 
defects. Conversely, from Eqs.(32) and (33) [and simi- 
larly from Eqs.(34) or (35)] it follows that, in the Uk — > 
limit, K n and V n do not tend to the values of the defect- 
free lattice [see, for instance, Fig. 1(b)]. This means that, 
in the non-Hermitian lattices synthesized in Sec. III. B, in- 
visibility of defects can be achieved. It should be noted 
that such non-Hermitian lattices with localized defects 
possessing a real- valued energy spectrum are not Vl~ in- 
variant, i.e. VT symmetry is not of relevance for the 
achievement of invisibility of the defects. 
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FIG. 5: (color online) Propagation of an initial Gaussian- 
shaped wave packet (snapshot of the site occupation proba- 
bilities P n (t) = |(7/)(f)|n}| 2 ) (a) in the non-Hermitian lattice 
with hopping rates shown in Fig. 1(b) [parameter values are 
given in the text], and (b) in the defect-free lattice. In (c) 
the behaviors of site occupation probabilities P n (t) at time 
t = 70 in the two lattices are depicted (the thin solid line 
refers to the non-Hermitian lattice with defects; the dashed 
line, almost overlapped with the solid one, refers to the defect- 
free lattice). The inset in (c) shows the behavior of the total 
occupation probability P(t) = 2~2 n Pn(t) versus time in the 
non-Hermitian lattice of (a). 



We have checked these predictions by direct numeri- 
cal simulations of wave packet propagation in Hermitian 
and non-Hermitian tight-binding lattices with zero site 
energies and with hopping amplitudes defined accord- 
ing to Eqs.(30) and (35), respectively. As an example, 
Fig. 3 (a) shows the propagation of an initial Gaussian- 
shaped wave packet \ip(t = 0)) = J2 n N exp{—(n + 
no) 2 /w 2 ] exp(— iqon)\n) in an Hermitian lattice with hop- 
ping rates given by Eq.(30) for parameter values k = 1, 
N = 3, U)\ = 0.6, a = 0, no — 70, wq = 10 and qo = n/2 
(J\f is the normalization constant). The profile of hop- 
ping rates for this lattice was shown in Fig. 1(a). For 
comparison, Fig. 3(b) shows the propagation of the same 
wave packet in the defect-free lattice. The distribution 
of site occupation probabilities P n (t) = \(n\ij](t))\ 2 at 
time t = 70 in the two cases is shown in Fig. 3(c). Note 
that, according to the previous analysis, the wave packet 
is fully transmitted in both lattices, and far from the 
inhomogeneities it propagates with the group velocity 
v g = 2Ksin(<7o) = 2. However, in the lattice with de- 
fects the wave packet is advanced, as one can see clearly 
from an inspection of Fig. 3(c). The behavior of the phase 
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(p(q) of the transmission coefficient of the partner lat- 
tice corresponding to the simulation of Fig. 3(a) is shown 
in Fig.f(c). One might think that, to make the Her- 
mitian lattice invisible, one should reduce the value of 
lu\; however, as discussed previously and as shown in 
Fig. 4, as uji is diminished toward zero, the defects in 
the hopping amplitudes vanish and the lattice basically 
becomes defect-free. Conversely, Fig. 5 shows that a non- 
Hcrmitian lattice can be invisible yet presenting defects 
in the hopping amplitudes. Figure 5(a) shows the prop- 
agation of the same initial Gaussian-shaped wave packet 
\tj)(t = 0)) = ^ n ATcxp[-(n + n ) 2 /w 2 }exp(~iq n)\n) 
of Fig. 3, but in the non-Hcrmitian lattice with hopping 
rates given by Eq.(35) for parameter values k = 1, N = 3, 
wi = 0.01 and a = 0.5 [the distribution of hopping rates 
for this lattice was shown in Fig. 1(b)]. For comparison, 
Fig. 5(b) shows the propagation of the same wave packet 
in the defect-free lattice. The distribution of site occu- 
pation probabilities P n (t) = \(n\ip(t))\ 2 at time t = 70 
in the two cases is shown in Fig. 5(c). Note that, ow- 
ing to the flatness of the phase ip(q) for this lattice [see 
Fig. 1(d)], the wave packet is fully transmitted with no ap- 
preciable delay and/or distortion, as one can infer from 
an inspection of Fig. 5(c). An outside observer thus can 
not distinguish whether the transmitted wave packet has 
been propagated in a defect-free or in an inhomogeneous 
lattice, and thus the defects in the non-Hermitian lat- 
tice are fully invisible. It should be finally noted that 
the total probability P{t) = J2 n P n {t) is transiently not 
conserved in the non-Hcrmitian lattice, and turns out to 
be amplified during interaction with defects, as shown in 
the inset of Fig. 5(c). Such an enhancement of the prob- 
ability, however, is not visible to the outside observer. 



V. CONCLUSIONS 



In this work we have investigated theoretically the is- 
sue of invisibility of rcflcctionless tight-binding lattices 
with defects synthesized by the intertwining operator 
technique of supersymmetric quantum mechanics. As 
for Hermitian lattices the defects are not invisible and 
time-of-fiight measurements of wave packets crossing the 
defect region may reveal their existence, in this work it 
has been shown that, in a certain class of non-Hermitian 
lattices with complex hopping amplitudes, the defects 
may appear fully invisible to an outside observer. In 
spite of non-Hermiticity, such lattices have a real-valued 
energy spectrum. As discussed in the Appendix C, 
arrays of evanescently-coupled optical waveguides with 
suitable longitudinal modulation of loss/gain coefficients 
and propagation constants could provide a physically re- 
alizable system to test invisibility in non-Hcrmitian tight- 
binding lattices. 
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Appendix A 

In this Appendix, the following theorem is proved: 
Let Hi be a tight binding Hamiltonian with = 0, 
and let H2N+1 a partner Hamiltonian synthesized from 
Hi by adding 2N new energy levels fj,i, ^2, M3, (J-2N, 
with [i2 = -fii, /i 4 = -fJ-3,-—, H2N = -H2N-1- Then 
H2N+1 can be constructed in such a way that Vn 2N+1 ^ = 
0. 

Let us first prove the theorem for N = 1. The partner 
Hamiltonian H2 is first constructed following the pro- 
cedure described in Sec. II, and the corresponding hop- 

(2) . . (2) 

ping amplitudes n n and site energies V„ are given by 

Eqs.(ll) and (12), respectively, with V^p =0. To syn- 
thesize the Hamiltonian H3, we need to construct the 

sequence </>^ satisfying the difference equation 



+ «. 



(2) A2) 
n+lVn+1 



+ V^^=^ (Al) 



and apply again the intertwining operator technique after 

the factorization H2 = Q2^-2+A*2- InEq.(Al), ^2 = — Mi 

(2) 

and the asymptotic behavior \<p n — > 00 for n — > ±00 
should be satisfied. A possible choice for the sequence 
cj>n can be obtained by observing that, since V™ = 0, 
from Eq.(2) it follows that |V>) = EnM)™^™ V) sat- 
isfies the equation Hi\ip) = —fXi\ip), and thus \<j)^) = 
1t\\ip) satisfies the equation %2\<i>^) = —fJ>i\((>^}, which 
is precisely Eq.(Al). Using Eqs.(6), (7) and (15) one 
obtains after some algebra 

0i 2) = -2(-l)"^W;. (A2) 

The hopping rates n n and site energies of the partner 
Hamiltonian W 3 = 7I2Q2 + A*2> obtained from % 2 after 
changing the order of the operators Q2 and IZ2 , are then 
given by [see Eqs.(ll) and (12)] 



(2) 

— ( 2 ) 



(A3) 



2 4 >(2 ' > (2) 

Ki 3 ' = ' + K n+1 , (2) K n ' s ( 2 ) (A4) 



x(2) 



where 



r (2) = _ 



A( 2 ) 



(A5) 
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Using in Eq. (A5) the expressions of defined by 

(2) 

Eqs.(6) and (11), and of </>„ as given by Eq.(A2), sub- 
stitution of Eq.(A5) into Eqs.(A3) and (A4) finally yields 
after some straightforward though lengthy algebra 



,(3) 



VP 



6 {1) 6 {1) 

Vn-lfn-2 



(A6) 
(A7) 



Therefore, for the partner Hamiltonian H3, obtained 
from Hi by adding the two energies [i\ and yU 2 = — Mi 



with the procedure described above, one has Vn 



(3) 



0. 



Starting from H3, one can repeat the procedure to con- 
struct a partner Hamiltonian H5 with = by adding 
to H3 the couple of eigenvalues /it 3 and = —[13. The 
hopping amplitudes n n of the new Hamiltonian will be 

given by Eq.(A6), with KrP and <jffl replaced by k^ 3 ' and 
(3) 

4>n , respectively. By induction, it follows that a part- 
ner Hamiltonian H 2 n+i, obtained from Hi by adding N 
couples of energies {^1,^2 = -fJ-i}, {^3,^4 = -^3}, 
{^2N-i, fJ-2N = —H2N-1}, can be always synthesized to 
have Vn 2N+1 ^ = 0, which proves the theorem. 



Appendix B 

In this Appendix we prove the Eqs.(21) and (22) given 
in the text relating the reflection and transmission coef- 
ficients of the two partner lattice Hamiltonians Hi and 
Hi- To this aim, let us first consider the case Hi > 2k, 



and let us indicate by lui the positive root of the equa- 
tion fii = 2KCOsh(wi). As — > k and — > at 
n — > ±00, the asymptotic behavior of 4>n\ satisfying 



Eq.(2), is of the form 



aexp(win) n — > +00 
/3exp(— win) n — > — 00 



where a and j3 are two non-vanishing constants. 
Eqs.(6), (7), (11) and (12) it then follows that 



.(i) 

n 

,(i) 

n 

,(2) 



— > — exp(=Fwi/2) for 
— > exp(±wi/2) for n 
— > 1 for n — > ±00 



->■ ±00 
±00 



K (2) -> for n-S- ±00. 



(Bl) 
From 

(B2) 
(B3) 
(B4) 
(B5) 



Let us then indicate by \ip^) = ^2 n V'™ \n) the solution 
to the equation Hi\tp^) = E^ 1 " 1 ) corresponding to the 
scattering of a forward propagating plane wave (coming 
from n — > — 00) with wave number q and energy E = 
2kcos((?) (0 < q < 7r). The eigenfunction IV'^) has 
therefore the asymptotic behavior expressed by Eq.(18) 
given in the text. According to Eq.(15), the function 

IV> (2) ) = E„tfV> with 



(2) 



,(lty(D 



(B6) 

satisfies the equation H2\ r 4>^) — E\ip^ 2 ' 1 ). Using 
Eqs.(18), (B2) and (B3), it follows that the asymptotic 
behavior of ipn^ is given by 



(2) 



exp(wi/2) + exp(— Wi/2 + iq)] cxp(-iqn) + ri{q) [exp(— oji/2 — iq) — cxp(wi/2)] cxp(iqn) n 
ti(q) [— cxp(c<Ji/2) + cxp(u;i/2 + iq)] cxp(-iqn) n 



00 



r 



_[B71 



i.e. \tp^) describes the scattering, in the lattice H2, of a 
plane wave with wave number q coming from n — > — 00 
and with amplitude [— exp(wi/2) + exp(— u>i/2 + iq)]. 
From Eq.(B7), the transmission (i 2 ) and reflection (r 2 ) 
coefficients of the partner lattice Hi are readily calcu- 
lated, obtaining the expressions (21) and (22) given in 
the text with 81 = 1. 

Let us now consider the case [ii < —2k, and let us 
indicate again by u>i the positive root of the equation 
2KCOsh(w 1 ) = — Hi. The asymptotic behavior of <j>n\ 
satisfying Eq.(2), is now of the form 



a(— 1)™ exp(win) n — > +00 
/3(— 1)" cxp(— coin) n — > — 00 



where a and (3 are again two non-vanishing constants. 
In this case, the asymptotic behavior of r^P and f„ \ as 
obtained from Eqs.(6), (7), and (B8), is given by 



n 



-iexp(=FWi/2) for n ->• ±00 (B9) 
-iexp(±wi/2) for rc ->• ±00. (BIO) 



As compared to the previous case \ii > 0, from Eqs.(B6), 
(B9) and (BIO) it follows that the asymptotic behavior 

/,(2) j 



(B8) of ip n is now given by the equation 
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(2) I — i [exp(wi/2) + cxp(— u>x/2 + iq)] cxp(— iqn) — ir\{q) [cxp(— uj\/2 — iq) + exp(o;i/2)] exp(iqn) n — > — oo 

-iti(q) [exp(— 0J\/2) + exp(o; 1 /2 + iq)} exp(-iqn) n — > oo 

fBll) 



r 



which replaces Eq.(B7). The transmission and reflection 
coefficients ti and r-i of the lattice H2 are readily calcu- 
lated from Eq.(Bll), and their expressions are given by 
Eqs.(21) and (22) with £1 = -1. 



Appendix C 

In this Appendix we briefly discuss a possible phys- 
ical realization of non-Hermitian tight-binding lattices 
with complex hopping rates, such as those discussed in 
Sees. III. B and IV. In the optical context, it is known 
that Hermitian lattices can be implemented by consid- 
ering light propagation in arrays of evanescently-coupled 
optical waveguides, the propagation direction z of light 
playing the role of time t in the quantum-mechanical 
problem (see, for instance, [Hi[26j]). The evolution along 
z of the modal amplitudes c n of light trapped in the var- 
ious waveguides of the array is governed by the tight- 
binding Hamiltonian (1), in which the site energies V n 
and hopping amplitudes be engineered by a suit- 

able design of waveguide channel widths, index changes of 
the guiding cores, and distances between adjacent waveg- 
uides in the array. In ordinary arrays, i.e. without loss 
or gain regions, V n and n n turn out to be real- valued, 
and thus the Hamiltonian % Hermitian. Non-Hermitian 
lattices with complex site energies can be mimicked by 
considering arrays of evanescently-coupled waveguides in 
which light propagation in each waveguide is either ab- 
sorbed or amplified by some loss or gain mechanism 
(see, for instance, [l9|, [25|), where the z-invariant gain 
or loss coefficients in the various waveguides determine 
the imaginary parts of the site energies V n . Such non- 
Hermitian lattices have been intensively investigated in 
the past few years, especially in connection with VT- 
symmetric quantum mechanics [TTMT9L [25jj . However, the 
non-Hermitian lattices that realize invisibility, discussed 
in Sees. III.B and IV, have real-valued site energies V n 
but imaginary hopping rates n n at some lattice sites. To 
implement in optics such invisible lattices, let us consider 
an array of evanescently-coupled waveguides and assume 
that a suitable longitudinal and periodic modulation of 
both gain/loss coefficient and effective modal index, with 
spatial period A, is impressed to some waveguides in the 
lattice. In this case, coupled-mode equations describing 
the evolution of the modal amplitudes c„ of light trapped 
in the various waveguides read (see, for instance, [26[) 

dc n 

i—p- = A„c„_i + A„ + ic n+ i + [V n + f3 n (z) - ij n (z)) c n 



where A„ is the (real-valued) coupling rate between 
waveguides n and n + 1, V n is the propagation constant 
mismatch from a reference value, and /3 n (z), j n (z) are 
the impressed longitudinal modulations of the propaga- 
tion constant and loss/gain coefficient, respectively. We 
assume that both /3 n (z) and "f n {z) are periodic functions, 
with spatial period A and with zero mean. This means 
that, on average, a light field propagating in a single 
waveguide of the array would not be damped nor am- 
plified. Assuming that the spatial period A of the mod- 
ulation is much shorter than the typical coupling lengths 
(~ 1/A„) and mismatch lengths (~ 1/V n ), after intro- 
duction of the amplitudes 

a n (z) = c n (z) exp [iip n (z)\ (C2) 

where 

y n {z)= f W»(£)-n»(0L (C3) 

Jo 

a set of effective equations for the slowly- varying ampli- 
tudes a n {z) can be derived by a multiple-scale analysis 
(see, for instance, They read explicitly 

.dan = A n (exp[i^„(z) -•i< / 9 n _i(z)])a„_i + (C4) 
dz 

+ A n+ i(exp[£^„(z) - iip n+ i(z)])a n+ i + V n a n , 

where (...) denotes the average with respect to z over the 
spatial oscillation period A. Let us then assume that: 

(i) /3 n (z) = p n P{z) and j n (z) = Pnl(z), where p n can 
take the values or 1. This means that some waveg- 
uides in the array are not modulated (those such that 
Pn = 0), whereas the modulated waveguides (those with 
p n = 1) have the same modulation profiles of loss/gain 
and propagation constant, defined by the two real- valued 
functions j(z) and P(z), respectively. 

(ii) The modulation functions j3(z) and 7(2) are chosen 
such that 

(exp[i<p(z)]) = (exp[-i<p(z)]) = iT, (C5) 

where <p{z) is defined by Eq.(C3) with f3 n = P(z) and 
7„ = "f(z), and T is a real- valued constant. 
Under such assumptions, Eqs.(C4) reduce to the follow- 
ing ones 

^dan _ Knan _ i _|_ Kn+ian+1 4. V„a n , (C6) 
dz 

where 

_ J A„ if Pn-l = p n 

Kn -\iTA n if Pn-X+Pn- 1 ' 
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and r is defined by Eq.(C5). In this way, Eqs.(C6) de- 
scribe the dynamics in a tight-binding lattice with hop- 
ping amplitudes K n between adjacent sites \n) and \n — 1) 
which can assume either real values (when the waveguides 
n and n— 1 are both modulated or both not modulated) or 
purely imaginary values (when one of the two waveguides 
n or n — 1 is modulated, but the other it is not). The ex- 
amples of reflectionless non-Hermitian lattices discussed 
in Sees. III. B and IV belong to such a class of lattices. It 
should be noted that satisfaction of Eq.(C5) requires a 
proper choice of the modulation amplitudes for loss/gain 
and propagation constant profiles. For instance, let us 
assume a sinusoidal modulation 

7 (z) = Ay cos(2ttz/A) , /3{z) = A fj cos(2ttz/A). (C8) 

In this case, from Eqs.(C3) and (C5) one obtains 



\ Z7T 



and hence the amplitudes Ap and A~ t must be chosen 
in such a way that the zero-order Bessel function Jo at 
the complex argument A(A ( g + iA g )/(2ir) gives a purely 
imaginary value. There are several possibilities to satisfy 
such a condition; for instance, one could fix the product 
ApA and determine, correspondingly, the product A y A-. 
for example, a choice can be 



2tt 



AA y 

~2~r 



2.096 



(CIO) 



(C9) which yields V ~ 1.941. 
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